Manifestation of triplet superconductivity in superconductor-ferromagnet structures 
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We study proximity effects in a multilayered superconductor/ferromagnet (S/F) structure with 
arbitrary relative directions of the magnetization M. If the magnetizations of different layers are 
collinear the superconducting condensate function induced in the F layers has only a singlet compo- 
nent and a triplet one with a zero projection of the total magnetic moment of the Cooper pairs on the 
M direction. In this case the condensate penetrates the F layers over a short length £j determined 
by the exchange energy J. If the magnetizations M are not collinear the triplet component has, in 
addition to the zero projection, the projections ±1. The latter component is even in the momentum, 
odd in the Matsubara frequency and penetrates the F layers over a long distance that increases with 
decreasing temperature and does not depend on J ( spin-orbit interaction limits this length). If 
5— I ■ the thickness of the F layers is much larger than £j, the Josephson coupling between neighboring 

S layers is provided only by the triplet component, so that a new type of superconductivity arises 
in the transverse direction of the structure. The Josephson critical current is positive (negative) for 
the case of a positive (negative) chirality of the vector M. We demonstrate that this type of the 
triplet condensate can be detected also by measuring the density of states in F/S/F structures. 
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Multilayered superconductor/ferromagnet (S/F) structures are under an intensive study now (for a recent review 
see e.g. 1 ). The interest in such systems originates from a possibility to find new physical phenomena as well from the 
hope to construct new devices based on these structures. Although a ferromagnet F attached to a superconductor S 
is expected to suppresses the order parameter in S, under certain conditions superconductivity and ferromagnetism 
may coexist and exhibit interesting phenomena. 

One of them is a nonmonotonic dependence of the critical temperature T c of the superconducting transition in S/F 
I/-) I multilayered structures on the thickness (If of the ferromagnetic layers. Theory of this effect has been developed in 
C<~) • Refs. 2 , and experimental results have been presented in Refs. 3 . 

Another interesting phenomenon is a 7r— state that can be realized in SFS Josephson junctions. It was shown 4 
that for some values of parameters (such as the temperature T, the thickness (If, the exchange energy J) the lowest 
Josephson energy corresponds not to the zero phase difference ip, but to ip = tt (negative Josephson critical current 
I c ). Detailed theoretical studies of this effect have been presented in many papers 7,8 . The 7r-state has been observed 
experimentally in Refs. 6 . 

Later it was discovered that the critical current I c in Josephson junctions with ferromagnetic layers is not neces- 
sarily suppressed by the exchange interaction and it may even be enhanced. Such an enhancement of I c has been 
demonstrated by the present authors on a simple model of a SF/I / FS junction, where / stands for a thin insulating 
layer 10 . It was shown for thin S and F layers that at low temperatures the critical current I c in a SF/I /FS junction 
may become even larger than in the absence of the exchange field (i.e. if the F layers are replaced by N layers, where 
J> , N is a nonmagnetic metal). More detailed calculations of I c (for arbitrary S/F interface transmittance) for this and 
similar junctions have been performed later in Refs. 11 . 

Properties of superconductors in S/F structures may change not only due to the proximity effect but also due 
■ to the long-range magnetic interaction. A spontaneous creation of vortices caused by the magnetic interaction has 
been predicted in a S/I/F system ( I is an insulating layer) 9 . In most papers on S/F structures the case of collinear 
(parallel or antiparallel) orientations of the magnetization M was considered. If the magnetization vector M is not 
constant in space, as in a domain wall, or if the orientations of M in different F layers are not collinear to each other, 
a qualitatively new and interesting effect occurs. For example, if a ferromagnetic wire is attached to a superconductor, 
a domain wall in the vicinity of the interface can generate a triplet component of the superconducting condensate 12 
(a similar case was analyzed in a later work 13 ). 

The existence of the triplet component (TC) has far reaching consequences. It is well known that the singlet 
component (SC) penetrates into a ferromagnet over the length £,/ = \J Dp / J , where Dp is the diffusion coefficient in 
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F. In contrast, it was shown that even for J >> T the TC penetrated F over a much longer distance £t = \/F>f/2ttT. 
This long-range penetration of the TC might lead to an increase of the conductance of the F wire if the temperature 
is lowered below T c 12,13 . 

In this paper we consider a multilayered S/F structure. Each F layer has a constant magnetization M but the 
direction of the M vector varies from layer to layer. We show that, in this case, the triplet component of the 
superconducting condensate is also generated and it penetrates the F layers over the long length £r that does not 
depend on the large exchange energy J at all. 

If the thickness of the F layers dp is much larger than £j, then the Josephson coupling between adjacent S layers 
and, therefore, superconductivity in the transverse direction is due to the TC. In the vicinity of the S/F interface 
the amplitudes of the SC and TC may be comparable but, unlike the TC, the SC survives in F only over the 
short distance £j from the S/F interface. In other words, in the multilayered F/S structures with a non-collincar 
magnetization orientation, a new type of superconductivity arises. The non-dissipative current within the layers is due 
to the s-wave singlet superconductivity, whereas the transversal supercurrent across the layers is due to the s-wave, 
triplet superconductivity. 

It is important to emphasize (see Ref. 12 ) that the TC in this case differs from the TC realized in the superfluid 
He 3 and, for example, in materials like Sr 2 Ru04 14 . The triplet-type superconducting condensate we predict here is 
symmetric in momentum and therefore is insensitive to non-magnetic impurities. It is odd in frequency and is called 
sometimes odd superconductivity. 

This type of the pairing has been proposed by Berezinskii in 1975 15 as a possible candidate for the mechanism 
of superfluidity in He 3 . However, it turned out that another type of pairing was realized in He 3 : triplet, odd in 
momentum p (sensitive to ordinary impurities) and even in the Matsubara frequencies w. Attempts to find conditions 
for the existence of the odd superconductivity were undertaken later in several papers in connection with the pairing 
mechanism in high T c superconductors 16 (in Ref. 16 a singlet pairing odd in frequency and in the momentum was 
considered). It is also important to note that while the symmetry of the order parameter A in Refs. 14-16 differs 
from that of the BCS order parameter, in our case A is nonzero only in the S layers and is of the BCS type. It is 
determined by the amplitude of the singlet component. Since the triplet and singlet components are connected which 
each other, the TC affects A in an indirect way. 

Therefore the type of superconductivity analyzed in our paper complements the three known types of super- 
conductivity: s-wave and d-wave singlet superconductivity that occur in ordinary superconductors and in high T c 
superconductors respectively, and the p-wave superconductivity with triplet pairing observed in Sr 2 Ru0 4 . 

In addition, the new type of the triplet superconductivity across the S/F layers shows another interesting property 
related to the chirality of the magnetization M. If the angle of the magnetization rotation 2a across the Sa layer 
(see Fig. 7) has the same sign as the angle of the M rotation across the Sb layer, then the critical Josephson current 
I c between Sa and Sb is positive. If these angles have different signs, then the critical current I c is negative and ir— 
state is realized (in this case spontaneous supercurrents arise in the structure). This negative Josephson coupling, 
which is caused by the TC and depends on chirality, differs from that analyzed in Ref. 4 . Depending on the chirality 
an "effective" condensate density in the direction perpendicular to the layers may be both positive and negative. We 
note that a dependence of the Josephson current on chirality has also been obtained in Ref. 5 . The authors of Ref. 5 
considered two magnetic superconductors S m with spiral magnetization, separated by a thin insulating layer I. In the 
latter case the TC exists in the bulk superconductors together with the SC ( and they cannot be separated) and the 
Josephson current depends on the chirality of the spiral structures. The main difference between our system and the 
system considered by the authors of Ref. 5 is that in our case only the long-range TC survives in the F layers whereas 
in the S m IS m junction both the SC and TC exist simultaneously. Therefore in the case of a collinear alignment of M, 
the Josephson coupling (and triplet superconductivity in the transverse direction) disappears in our system, whereas 
it remains in the S m IS m system. 

Another possible detection of the TC in the S/F structures may be achieved by measuring the density of states 
(DoS) in a F/S/F trilayer (see Fig. 1). We will see that the long-range TC causes a measurable change of the local 
DoS at the outer side of the F layers even if (If is much larger than £j. 

The plan of this paper is as follows. In the next section we make some preliminary remarks concerning the TC in 
S/F structures. We consider a three-layer FSF structure and calculate the condensate function in this structure. We 
show that the amplitude of the TC is proportional to sin a and its long-range part is an odd function of the Matsubara 
frequency w (the SC is an even function of w), where ±a is the angle between the z-axis and the magnetization in 
the right (left) F layers. We discuss properties of the TC and calculate the DOS related to it. In Sec. 3 we calculate 
the Josephson current between adjacent 5* layers and discuss its dependence on the chirality of the magnetization 
variation in the system. In Sec. 4 we take into account spin-orbit interactions and study the effect of this interaction 
on the TC. In the conclusion we discuss the obtained results and possibilities of an experimental observation of the 
predicted effects. The odd triplet superconductivity in F/S structures was first predicted by the present authors in 
a short paper where the case of small angles a and of a perfect F/S interface was considered 17 . 
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II. THE CONDENSATE FUNCTION IN A F/S/F SANDWICH 



In order to get a better understanding of the properties of the superconducting condensate in the presence of the fer- 
romagnetic layers, we consider in this section a simple case of a trilayered F/S/F structure (see Fig.l). Generalization 
to a multilayered structures is of no difficulties and will be done in the next section. 

In the most general case, when the magnetization vectors M of the F-layers are non-collinear, the electron Green 
functions are 4x4 matrices in the particle-hole<g)spin-space 18 . The 4x4 matrix Green functions have been introduced 
long ago 20 and used in other papers 21 . Later on they were used in Ref. 22 for a description of magnetic superconductors 
with a rotating magnetization. 

A very convenient way for the study of proximity effects is the method of quasiclassical Green's functions 23-25 . 
Equations for the quasiclassical Green's functions have been generalized recently to the case of a non-homogeneous 
exchange field (magnetization) M 26 . 

Following the notation of Ref. 8 we represent the quasiclassical Green functions in the form 
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^£fo)»n" I ^(^« S (t)i + V (0) , (1) 



where the subscripts n and s stand for the elements in the particle-hole and spin space, respectively, and T3 is the 
Pauli matrix. The field operators tp ns are defined as tjj\ s = ip s and ip2s — tpf ( s denotes the opposite to s spin 
direction) . 

The elements of the matrix g diagonal in the particle-hole space (i.c proportional to f and T3) are related to the 
normal Green's function, while the off-diagonal elements (proportional to T\ and T2) determine the superconducting 
condensate function /. In the case under consideration the matrix (1) can be expanded in the Pauli matrices in the 
particle- hole space (fo is the unit matrix): 

9 = %t q + .9 3 ^3 + / , (2) 

where the condensate function is given by 

/ = Am + f 2 ir 2 . (3) 

The functions g i and f i are matrices in the spin-space. In the case under consideration the matrices j i can be 
represented in the form 

h(x) = fo(x)a + f 3 (x)& 3 (4) 

f 1 (x)=f 1 (x)a 1 (5) 

This follows from the equation that determines the Green's function (see below). 

Let us discuss briefly properties of the condensate matrix function /. According to the definitions of the Green's 
functions, Eq. (1), the functions fi(x) are related to following correlation functions 

h ~ (V'tV'i) - 0>li>l) , 

f ~ (V T Vi) + (V>iV>T> , (6) 



-a 



a 



-(d s + d F ) -d s d s d s + d F x 

FIG. 1. The F/S/F trilayer. The magnetizations vectors in the F layers make an angle ±a with the z-axis, respectively 
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The function / 3 describes the SC, while the functions /o and fi describe the TC (see for example Ref. 27 ). The 
function /g is proportional to the zero projection of the triplet magnetic moment of the Cooper pairs on the z-axis, 
whereas the function fi corresponds to the projections ±1. 

It is important that in the absence of an exchange field J (or magnetization M) acting on spins, the SC, i.e. the 
function / 3 , exists both in the superconducting and normal (non-magnetic) layers. If J is not equal to zero but is 
uniform in space and directed along the z-axis, then the part /o of the TC arises in the structure. 

However, both the functions / 3 and /o decay very fast in the ferromagnet (over the length The singlet 

component decays because a strong magnetization makes the spins of a pair be parallel to each other, thus destroying 
the condensate. The triplet component with the zero projection of the magnetic moment is also destroyed because it 
is more energetically favorable for the magnetic moment to be parallel to the magnetization. 

On the other hand, the structure of the matrix / (the functions /J depends on the choice of the z-axis. If the 
uniform magnetization M is directed not along the z-axis (but, say, along the x-axis), terms like f^iri inevitably 
appear in the condensate function (see for example, Ref. 22 where such a term was obtained even at Q = 0, Q is 
the wave vector of a spiral magnetic structure). However, the condensate component corresponding to this term 
penetrates the F layer over the short distance £j only. 

Therefore, we can conclude that the presence of terms like f y%T\ in the condensate function does not necessarily 
mean that the TC penetrates the F layer over the long distance £t- Actually, long-range effects arise only if the 
direction of the vector M varies in space. If the magnetization has different directions in neighboring F layers, then 
not only /o but also /i arise in the system and both functions penetrate the ferromagnetic layer over a long distance 

In order to find the Green's function g, we consider the diffusive case when the Usadel equation is applicable. This 
equation can be used provided the condition Jt <C 1 is satisfied ( r is the momentum relaxation time). Of course, 
this condition can hardly be satisfied for strong ferromagnets like Fe and in this case one should use a more general 
Eilenberger equation for a quantitative computation. However, the Usadel equation may give qualitatively reasonable 
results even in this case. 

The Usadel equation is a nonlinear equation for the 4x4 matrix Green's function g and can be written as 

Dd x (gd x g) - uj [f 3 <7 , <?] + iJ {[f 3 <7 3 ,<?] cosa(x) + [f <T 2 ,g] sino:(x)} = -i [A,g] . (7) 

In the S layer D = D$, J — 0, A = AiT 2 <T 3 (the phase of A is chosen to be zero). In the F layers D = Dp, a(x) = ±a 
for the right (left) layer and A = 0. Eq. (7) is complemented by the boundary conditions at the S/F interface 28 

l(g9 x g) F = (gd x g) s , x = ±d s (8) 
27bO (gdx9) F = ± [g s , 9f\> x = ±d s , (9) 

where 7 = apjas, and <js,f are the conductivities of the F and S layers, 7& = crpRb/O is a coefficient characterizing 
the transmittance of the S/F interface with resistance per unit area Rb- 

If linearized, the Usadel equation can be solved analytically rather easily. The linearization may be justified in the 
two limiting cases: a) T is close to the critical temperature of the structures T* (the latter can be different from the 
critical temperature of the bulk superconductor T c ), and b) the resistance of the S/F interface Rb is not small. In the 
latter case the condensate function in the S layer is weakly disturbed by the F film and the function / 3 in Eq. (3) 
can be represented in the form 

h(x) = fs + Sf 3 (x), \x\<d s , (10) 

where fs = A/iE^ and E u = ^lu 2 + A 2 . The function <5/ 3 as well as the functions fo,i are assumed to be small. In 
the F layers all the components of the condensate function / are small. The functions g and g 3 in Eq. (2) in the 
superconductor are given by 

<?3 = ^o(.gs + <5.9o) + 0-353 (11) 
9o = °2.92 • (12) 

Here gs = sgnoj.gs, gs = M/E^. From the normalization condition 

f = 1 (13) 
we obtain expressions relating the functions 5go, 32,3 to the functions <5/ 3 , /o,i 

Sgo = {fs/gs)Sf 3 , 93 = (fs/gs)fo, 92 = {fs/gs)h ■ (14) 
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Now we linearize Eq. (7) with respect to Sf = iT2(a 3 5f 3 + &ofo) + iTi&ifi an d obtain 

d 2 x jf-4sf = o 

in the S layer, and 

d lx 5 f ~ K l S f + iK 2 j{T [a 3 , Sf} + cosa ± f 3 [<7 2 , Sf]- sina} = 



(15) 



(16) 



in the F layers. Here k% = 2E U /D S , «£ = 2\w\/D F , Kj = Jsgncu/Dp and [A,B}± = AB ± BA. The signs ± in Eq. 
(16) correspond to the right and left layer respectively The corresponding linearized boundary conditions for Sf arc 



(l9s)d x f F = d x Sf s 

±lbtjdf F = -(fs + Sf s ) + g s f F , 



(17) 
(18) 



where f s — iT 2 a 3 fs and the signs ± correspond to the right and left layer. Solutions for Eqs. (15-16) can be written 
as a sum of exponential functions exp(±Kx), where the k's are the eigenvalues of Eqs. (15-16). In the S layer the 
equations for Sf 3 , /o,i are decoupled and there is only one eigenvalue k = ks- In the F-layers the equations are coupled 
and there are three different eigenvalues 17 



Kl,2 = k± ~£/(l±i) , 
K3 = K w = y/2\uj\/D F . 



(19) 
(20) 



We see from these equations that two completely different lengths £j and £t determine the decay of the condensate 
in the F layers. At all temperatures T < T* the length £x much exceeds £j and is the same the length describing the 
decay of the standard singlet condensate in a normal metal. 

We have assumed that J » T*, which is realistic unless the exchange field is extremely small. In order to find 
analytical expressions for the functions fi we also assume that the thicknesses of the S and F layers satisfy the 
conditions 



d S « £s = VD S /2ttT*, d F »£j. 

In this case the solutions for Eqs. (15-16) have the form 

Sh{x) = a 3 cosh(K S x) 
fo(x) = a cosh(K S x) 
fi(x) = ai sinh(K S x) , 



in the S layer and 



coshK w (x - d s - d F ) 
fi (x) — bi — — h sgnw sin a 



cosh(K^dp) 

coshK w (a; — ds — dp) 



fo(x) = - tana hi 

cosh^udF) 

h(x) = & 3+e - K +(*- d =) + b 3 -e- K -( x - 



+ sgnw cos a 



-b 3+ e- K+{x - da) +b 3 _e- K - {x - d ° ) 



s) 



(21) 



(22) 
(23) 
(24) 



(25) 

(26) 
(27) 



in the right F layer. The solutions in the left F layer can be easily obtained recalling that the function fi(x) is odd 
and fo,3(x) are even functions of x. From Eqs. (22-27) and the boundary conditions Eqs. (17-18) we find 



fs 



ks tanh @sM^ 



M+T- 



-M_T+ 

k s (k + — K_)sgnaj 



h = M 5 s + 7^tanh6.) = -fs^ cosh , es{M+T +M T+) , 

where 6 S = K s d s , ®f = ^dp, k± = K±/(g s + 7fcO«±), « = K u /(g s + 7&£/«u> tanhG^), k s = n s /(gsl) and 

M± = T± (ks coth 65 + k tanh Op) + tan 2 a C± (ks tanh 6s + k tanh 6 p) 
T± = ks tanh 6 s + «± 
C± — ks coth 6s + k± . 



(28) 
(29) 
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The solutions presented above are valid if the correction 5/3 to the condensate function f s in the S layer is small 
(in the F layer Sf 3 is even smaller ). From Eqs. (17-18) one can readily see that the condition 

5fa(d s ) ~ 5f 3 (0) = a 3 cosh9 s = b 3+ + b 3 _ -f s <l (30) 

should be satished. Here |0g| <C 1 is implied. Actually we have neglected the term 6f 3 in the normalization condition 
(13) assuming that 5f 3 <C 1 (see Fig. 2). 

The amplitude 03 of the SC depends on many parameters, such as temperature (energy), 7^, etc. Therefore, the 
validity of our approach should be checked for every set of parameters. If we are interested in thermodynamical 
quantities such as the critical temperature or the Josephson current, we may set tu ~ max{T, A}. When calculating 
the density of states the situation is different because /s(e) has a singularity at e = A which is rounded off by a 
damping factor in the quasiparticle spectrum. In this case our approach breaks down near the energy e ~ A (see 
Fig. 2), when the condition (30) is violated. It is also clear that our approach is valid provided either the temperature 
is close to the critical temperature T* of the system or 7f, is not too small. 




0.5 1 1.5 2 2.5 
8/A 

FIG. 2. Dependence of \a§\ 2 on the energy e. The dashed vertical lines show the region in which our approach fails. Here 
7 = 0.05, J/A = 25, ds/^A = 0.4, cZf/£a = 0.5, 7b = 0.5, a — 7r/4 and the damping factor F — 0.1. We have defined 
£a = \J Ds /A, where A is the BCS order parameter. 

Now we discuss the properties of the obtained solutions (Eqs. (22-29)). From Eqs. (27-28) one can see that the SC 
is an even function of to and decays sharply in the ferromagnet over the short distance £j. In contrast, the amplitudes 
of the TC fo and f\ are odd functions of u> and penetrate the ferromagnet over the longer distance £t = \/Dp /2-kT. 
The long-range part of TC determined by the amplitude b\ has the maximum at a = 7r/4. This value of a corresponds 
to a perpendicular orientation of the magnetizations in the F layers. For a parallel (a = 0) or antiparallel alignment 
of the magnetizations (a = n/2) this amplitude decays to zero. In Fig. 3 we plot the spatial dependence of the SC 
and the long-range part of the TC. We see that both amplitudes are comparable at the S/F interface but the SC 
decays faster than the TC. 




-(ds+dF) -ds ds ds+dF 

FIG. 3. The spatial dependence of Im(SC) (dashed line) and the long-range part of Re(TC) (solid line). We have chosen 
7 = 0.2, J/T c = 50, 76 = 0.05, dpy/Tc/Ds = 2, d,s sjTcjDs = 0.4 and a = vr/4. The discontinuity of the TC at the S/F 
interface is because the short-range part is not shown in this figure. 
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The long-range part of TC leads to interesting observable effects that will be discussed in the next sections. In 
Refs. 12,13 the conductance of a ferromagnetic wire attached to a superconductor was calculated. It was assumed that 
the F wire had a domain wall located at the S/F interface. This inhomogencity of the magnetization induces a TC, 
which leads to an increase of the conductance for temperatures below T c . 



A. Critical temperature 



In this section we discuss briefly the effect of the TC on the critical temperature T* of the structure. For the 
parallel and antiparallel alignment of the magnetizations the critical temperature of the multilayered structure T* 
was calculated in many papers 2,29 . The angle dependence of the critical temperature in a F/S/F structure was 
analyzed in Ref. 30 . However the form of the condensate function presented in Ref. 30 is not correct because the 
authors started from an equation different from Eq. (7). As a result, the long-range TC was completely lost. 

The equation that determines T* has the form (we assume that ds -C £s, see Refs. 2 ) 

v L ' w= — oo ^ ) 

We have obtained a solution for b$±, (Eq. (28)), assuming that A is constant in space ( this approximation corresponds 
to the so-called single-mode approximation used in many earlier works 2 ). It is established in Ref. 29 that for some 
parameters this approximation gives a rough estimate for T*. A careful analysis of Ref. 29 shows that T* remains 
finite even for values of the parameters 7, 76, kj, for which other approaches predict a zero critical temperature. We 
will not discuss quantitatively the dependence of T* on the angle a. Note however that, as follows from Eqs. (28) and 
(31), the critical temperature T* depends on a and cIf even in the case when dp 3> £j (if a 7^ 0). This dependence 
is due to the long-range part of the TC and, in order to determine it, one has, generally speaking, to go beyond the 
single-mode approximation. Note, however, that this dependence may be weak. 



B. Local density of states 



In this section we calculate the change of the local DoS in the F layers due to the TC. It is clear that, for distances 
from the S/F interface larger than £j, only the TC leads to a variation of the local DoS. Thus, if the thickness dp is 
much larger than £j one can detect directly the presence of the TC performing measurements of the DoS at the outer 
side of one of the F layers. Any deviation from the normal value would be only due to the TC. 

We calculate the local DoS at x = ds + dp. The expression for the normalized DoS is (we ignore the difference in 
the DoS for the up and down spin directions. This approximation is consistent with the quasiclassical assumption 
that J <C ep, where ep is the Fermi energy) 

u=-= l -Tr(f,a ){g R -g A ) , (32) 

Vq 8 

where v is the DoS in the normal state, thus v = 1 + 5v (Su is a correction due to the proximity effect). As it was 
mentioned before, in the case dp ^> £j only the TC (i.e. the functions fo(x) and fi(x)) contributes to the DoS. From 
the normalization condition Eq. (13) and Eq. (32) we obtain 

1 (b R ) 2 

5v = -Re — - v u - , (33) 

2 cos 2 a cosh Of 

where Qp — \J '— 2ie/ 'Dpdp, and b R is the amplitude of the retarded Green's function in Eqs. (25-26). It is obtained 
from b\ by replacing ui by — ie. In Figs. 4, 5 and 6 we plot the dependence of 5v on e for different a, c?f and 7;,, 
respectively. For the range of parameters chosen in these plots the function |as(e)| 2 has the shape shown in Fig. 2. 
Thus, our approach is valid almost for all energies and fails only in a very narrow region close to e = A. In order to 
avoid singularities in f R we have taken into account a finite damping factor T = 0.1 in the expression for f R : 
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As follows from Eq. (29) 8v is zero for a — 0, 7r/2. The largest change in the DoS is achieved when a = 7r/4 
(perpendicular orientation of magnetizations in the F layers). We see that the correction to the DoS is small but 
observable. Kontos et al. presented in Ref. 31 measurements of 8v in thin F layers (few nanometers). The order of 
magnitude of the observed 8v (~ I0~ 3 ) is the same as the presented in Figs. 4-6. However, in Ref. 31 the variation 
of the DoS was caused by the penetration of the SC into the F l ayer over t he short distance £j. In our case such a 
variation can be observed in much thicker F layers (dp ~ £t = ^/Df/^ttT* 3> £j)- 

It is interesting to compare our result for the FSF structure with non-collinear magnetization with corresponding 
results for NSN structures (N is a normal layer). At first glance, the behavior of the odd triplet condensate in the 
ferromagnet is very similar to that of the conventional singlet condensate in a normal metal. In both cases the 
amplitude of the condensate decays exponentially with the length £x ( Eq. (20)). However, there is an essential 
difference. In the N layer an energy gap is induced due to the singlet condensate. The value of the energy gap is 
determined by mm{ A, Djy / (aw cIn Rb)} 32 ■ In contrast, no subgap appears in the ferromagnet due to the triplet odd 
condensate considered here, although the TC penetrates over the entire F layer provided its thickness dp is not very 
large, dp < £r. The main reason for the absence of a subgap e sg in the FSF system is the following. In SN structures 
the condensate function is not small at energies |e| < e sg . The exchange field shifts this energy interval by the large 
value J so that at low energies the condensate function ( both singlet and triplet) is small if 75 is not too small. Note 
also that the amplitude of the TC is smaller than the amplitude of the SC in a NSN structure since it contains a large 
parameter k± ~ yj in the denominator (see Eqs. (28-29)). 

For completeness we finally note that the change of the local DoS in the ballistic case (Jt ^> 1) was considered in 
Ref. 33 and in the pure ballistic case (t — * 00) in Ref. 34 . It turns out that the results in these two cases differ greatly 
from those obtained in the present paper for a diffusive system (Jt <C 1). 

5v . . . . . , 

0.002 







-0.002 




FIG. 4. The normalized DoS 8u as a function of the energy for a = 37r/8 (solid line), a — n/8 (dashed line) and a = n/4 
(point-dashed line). Note that for a = 0,7r/2 8v = 0. We have chosen 7 = 0.05, J/A = 25, 75 = 0.5, g!f/£a = 0.5, and 
rfs/^A = 0.4. Here £a = yDs/A and A is the BCS order parameter. 
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0.002 




FIG. 5. The normalized DoS &v as a function of the energy for c(f/£a = 0.8 (solid line), c£f/£a = 1.2 (dashed line). The 
point-dashed line shows the contribution to the DoS from the SC (/a). The latter is multiplied by a factor of 100. We have 
chosen a — 7r/4. All other parameters are the same as in Fig. 4 



Sv 




FIG. 6. The normalized DoS Sv as a function of the energy for 75 = 0.5 (solid line), 75 = 1 (dashed line) and 75 = 1.5 
(point-dashed line). We have chosen (1f/(,a = 0.5. All other parameters are the same as in Fig. 5 

III. JOSEPHSON CURRENT IN A F/S/F/S/F STRUCTURE 

In this section we calculate the Josephson current between the S layers of a FSFSF structure. We assume again 
that the thickness of the F layers <If is much larger than £,/ (Eq. (21)). In this case the Josephson coupling between 
the S layers is due to the long range part of the TC. Therefore the supercurrent in the transverse direction is unusual, 
since it is caused by the triplet component of the condensate that is odd in frequency and even in momentum. 

At the same time, the in-plane superconductivity is caused mainly by the ordinary singlet component. Therefore 
the macroscopic superconductivity due to the Josephson coupling between the layers is an interesting combination of 
the singlet superconductivity within the layers and the odd triplet superconductivity in the transversal direction. 

We will see that the unusual character of the superconductivity in the transversal direction leads to peculiarities of 
the Josephson effect. For example, if the bias current flows through the terminal superconducting layer So and Sa 
(see Fig. 7), the supercurrent is zero because of the different symmetry of the condensate in So and Sa- In order to 
observe the Josephson effect in this structure the bias current has to pass through the layers Sa and Sb, as shown 
in Fig. 7. The supercurrent between Sa and Sb is non-zero because each superconductor has its "own" TC and the 
phase difference if is finite. 







So 



-a t a 



>B 



3a (pos. chirality) 
—a (neg. chirality) 



FIG. 7. The multilayered structure considered. The arrows show the bias current. In the case of positive (negative) chirality 
the magnetization vector M of the layer F3 makes an angle 3a (—a) with the z- axis, i. e. in the case of positive chirality the 
vector M rotates in one direction if we go over from one F layer to another whereas it oscillates in space in the case of negative 
chirality. 



The Josephson current Is is given by the expression 

Is = (L y L z )a F Tr (f 3 <7 ) £ faf ( 35 ) 

This current was calculated for the case of small angles a in Ref. 17 . Here L y L z is the area of the interface and of 
is the conductivity of the F layer. The simplest way to calculate the Is is to assume a weak coupling between the S 
layers, which corresponds to the case when the condition d F > £,t holds. In this case the long-range part of the TC 
is given by the sum of two terms each of those is induced by the layers Sa and Sb in Fig. 7: 

/» = f A (x) + S.U.f B (x -d s - d F )U+.S+ , (36) 

where 

/ A (a;) = e""-^-^) (h.ihti + b .iT 2 o ) (37) 

is the long range part of the TC induced by the layer Sa- The coefficient 61 is given by Eq.(29) and bo = — tan a 61. 
If the Sa,b/F interfaces are identical as well as the superconductors Sa and Sb , the function / B is equal to / A if one 
replace the exponential function exp(— k^(x — dg)) by exp(/-c w (a; — d$ — dp))- The phase of the Sa layer is set to be 
zero and the phase of the Sb is f. This phase has been taken into account by the gauge transformation performed 
with the help of the matrix S — f cos(</?/2) + if 3 sin(y>/2). The magnetizations M of the layers Fj and F 2 make an 
angle =Fa with the z-axis respectively. For the direction of M in the F 3 we consider two cases: a) the direction of 
magnetization is —a (negative chirality) or b) 2a (positive chirality). In the latter case the matrix U in Eq. (36) is 
given by 

U = T0&3 cos a + iT3&2 sin a . (38) 

In the case of negative chirality, U is the unit matrix and one has to change the sign of a in the expression for the 
function f B (Eq. (36)). In Fig. 8 we show schematically the spatial dependence of fi(x). 
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FIG. 8. The spatial dependence of the amplitude of the TC fi(x) in the case of positive (solid line)and negative (dashed 
line) chirality. 

Substituting Eq. (36) into Eq. (35) one obtains after simple transformations Ig = I c sin ip, where 

eR F I c = ±2TtT^K UJ d F b\{a) (l + tan 2 a) e ^ dF ^ , (39) 

a? 

where bi(a) is given in Eq. (29) and the sign "+" ("-") corresponds to the positive (negative) chirality. In the case 
of negative chirality the critical current is negative (7r-contact). It is important to emphasize that the nature of the 
7r-contact differs from that predicted in Refs. 4 and observed in Ref. 6 . In our case the negative Josephson coupling is 
due to the TC and can be realized in S/F structures with negative chirality. This gives a unique opportunity to switch 
experimentally between the and 7r-contacts by changing the angles of the mutual magnetization of the layers. It is 
worth mentioning that another effect concerning the chirality of the M vector was studied by the authors in Ref. 35 . 
It was shown that the resistance of a multi-domain ferromagnetic wire depends on the chirality of the M variation in 
space. 

In Fig. 9 we plot the dependence of I c on the angle a. If the orientation of M is parallel (a = 0) or antiparallel 
(a = 7r/2) the amplitude of the triplet component is zero and therefore there is no coupling between the neighboring 
S layers, i.e. I c = 0. For any other angle between the magnetizations the amplitude of the TC is finite. This leads to 
a non-zero critical current. At a = ir/A ( perpendicular orientation of M) I c reaches its maximum value. 




7c/4 71/2 

a 



FIG. 9. Dependence of the critical current ( normalized with respect to the maximum value) on the angle a. We have chosen 
the same values as in Fig. 4 

The weak coupling assumption (dp > £y) leads to an exponential decay of I c with increasing dp ( Eq. (39)). In 
the case dp < £r, Eq. (39) is not valid. One can easily obtain I c for the case of an arbitrary d F and small a. It turns 
out that in this case Eq. (39) remains valid if the exponential factor exp(— K u dp) is replaced by cosh -2 (n u dp / '2) and 
in the expression for b\ (Eq. (29)) Q F is replaced by Q F /2. 

In order to estimate the value of the critical current I c , we use Eq. (39). If dp exceeds the length £t ( for example 
d F I £t = 2) only the term with n = (i.e. u> = 7rT*) is important in the sum. In this case one obtains 

T* it \T*J V ' 

where the factor C can be easily expressed in terms of M±, T±, etc. Thus, C depends on many parameters such as 
7, 7^, Kj, etc. We estimate C for values of these parameters similar to those which were used in Ref. 29 : 7b = 0.5, 
7 = 0.1, d s n s = 0.4, dpK^ = 1.5, k^/ks = 3. We get C = 10~ 2 - 10~ 3 for njd s = 5-10. The expression (40) for I c 
also contains the parameters (A/T*) 2 and exp(— dpKr) which are also small. We note however that if dp < £t, the 
exponential function is replaced by a numerical factor of the order of 1. The factor (A/T*) 2 is also of the order 1 if the 
temperature is not close to T*. Taking Op 1 =60/if2.cm (c/. Ref. 29 ) and dp ~ £t ^200 nm we obtain I c ~ 10 4 — 10 5 
A. cm -2 ; that is, the critical current is a measurable quantity ( see experimental works 6 ) and the detection of the TC 
is possible. 



11 



IV. EFFECT OF SPIN-ORBIT INTERACTION 



So far the only interaction we have considered in the ferromagnet is the exchange field J acting on the conducting 
electrons. However, in reality spin-orbit interactions that appear due to interactions of electron spins with spin 
orbital impurities may become important. Following again the notation of Ref. 8 we write an additional term in the 
Hamiltonian which describes the spin-orbit part as 36 ' 37 

H '° = & E c ns P (p x P ') {s):?; Cn , s , pl , (4i) 

n.s,p,n ,s ,p 

where S = (<7i, 02, f 3 a 3 ) and p and p' are the momenta before and after scattering at the impurities. Although in 
general the characteristic energy of the spin-orbit interaction is much smaller than the exchange energy, it can be 
comparable with the superconducting gap A and therefore this effect should be taken into account when describing 
the supercurrent. 

In the Born approximation the self-energy is given by 

S so = n \U so \ 2 < G > s .o. , where < G > s . = v j d^ p J ^- (n x n') SGS (n x n') . (42) 

Here n is a unit vector parallel to the momentum. Including this term in the quasiclassical equations is straightforward 
and the resulting Usadel equation takes the form 36 

-iDd r (gd T g) + i (f 3 d t g + dfjfo) + [A, g] + J [n, g] + — [Sf 3 gf 3 S, g] = , (43) 

Ts.o. 

where 

— = \imn [ ^\U so \ 2 sm 2 9 (44) 

is the spin-orbit scattering time. 

As before, one can linearize Eq. (43) in the F-layer and obtain equations for the condensate function / similar to 
Eqs. (15-16) but now including the spin-orbit interaction term. The solution again has the form 

f(x) = if 2 ® (/o(z)<xo + j*3 (2)03) + ir\ ® fi{x)&i . (45) 
The functions fi(x) are given by fi(x) = bj exp[njx], where the new eigenvalues Kj are 

(46) 
(47) 

We see from these equations that the singlet and triplet components are affected by the spin-orbit interaction making 
the decay of the condensate in the ferromagnet faster. In the limiting case 4/r so > J, T c both the components 
penetrate over the same distance £ S . D . = \Jt so Df and therefore the long-range effect is suppressed. In this case the 
characteristic oscillations of the singlet component are destroyed 37 . In the more interesting case 4:/t so ~ T c < J, 
the singlet component is not affected and penetrates over distances of the order £j. At the same time, the triplet 
component is more sensitive to the spin-orbit interaction and the penetration length equals min(^ so ,^T)> £j. 

Spin-orbit interaction is relevant in systems with large Z elements. The characteristic spin-orbit energy l/r s . . also 
depends on scattering concentration and density of states (cf. Eq. (44)). Experimental dataconcerning this energy is 
still unclear and controversial, mainly due to the difficulty to separate the contribution of the spin-orbit from other 
scattering types. From numerical band structure calculations one can estimate the parameter Jt s o . For example, for 
a typical magnetic transition metal, like Fe, in the dirty limit Jt s . <~ 10 2 , while for dirty Gd Jt s . <~ 10 ( see Ref. 38 
and references therein). Thus, according to our model, material like transition metals are better candidates in order 
to observe the predicted effects. 

Thus, provided the spin-orbit interaction is not very strong, the penetration of triplet condensate over the long 
distances discussed in the preceding sections is still possible, although the penetration length is reduced. 
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V. CONCLUSION 



We studied odd, s-wave, triplet superconductivity that may arise in S/F multilayered structures with a non-collincar 
orientation of magnetizations. 

It was assumed that the orientation of the magnetization is not affected by the superconductivity (e.g. the energy 
of the magnetic anisotropy is much larger than the superconducting energy). The analysis was carried out in the dirty 
limit (Jt <C 1) when the Usadel equation is applicable. 

It was shown that for all values of a the condensate function consists of a singlet (SC) and a triplet (TC) components. 
Even in the case of a homogenous magnetization (a = 0), in addition to the SC, the TC with the zero projection 
onto the z axi s arises. In this case, both the SC and the TC decay in the F layers over a short distance given by 
£j = y/Dp/J. If the magnetization vectors M are not collinear a ^ 0,7r/2, all projections of the TC appear, in 
particular, those with non-zero projection on the z-axis. In this case, the TC penetrates the F layer over a long 
distance £t = \/ Dp/2irT. In the presence of spin-orbit interaction this penetration length is given by min(£ so ,£r), 
where £ so = \/T so Dp. Generally, this length may be much larger than £j. 

Thus, if the condition cIf 3> £j is fulfilled the Josephson coupling between neighboring S layers is only due to the 
TC. Therefore in this case a new type of superconductivity may arise in the multilayered structures with non-collincar 
magnetizations. The supercurrent within each S layer is caused by the SC, whereas the supercurrent across the layers 
is caused by the triplet condensate, which is odd in the frequency ui and even in the momentum. 

The TC in our case is completely different from the triplet condensate found in Sr 2 Ru04 14 . In the latter case one 
has a p-wave, even in ui, triplet superconductivity, which is suppressed by impurity scattering. In contrast, the TC 
we have considered is not affected by non-magnetic impurities. The reason for the existence of the long-range TC is 
the fact that if a 0, the SC and the TC are coupled and, in addition to k± — ± i), the eigenvalue Kp — C^ 1 

appears. The latter corresponds to the long-range penetration of the TC in the ferromagnct. 

The triplet superconductivity in S/F structures possesses an interesting property: the Josephson current depends 
on the chirality of the magnetization M: If the M vector rotates in only one direction (the positive chirality) the 
critical current I c is positive. If the direction of the M vector oscillates in space (the negative chirality) then I c < 0. In 
the latter case spontaneously circulating currents must arise in the structure. This result can be explained as follows: 
if the chirality is positive the averaged M vector < M > is zero and the S/F structure behaves as a superconductor 
with anisotropic properties (the singlet superconductivity along the layers and the triplet superconductivity across 
them). In the case of the negative chirality the average in space yields a non-zero magnetization < M >^ 0. In such 
a superconductor with a build-in magnetic moment the circulating currents arise as they arise in superconductors of 
the second type in the mixed state. 

Note also that in a single Josephson FSFSF junction a non-zero magnetic field exists also inside the junction and 
this causes Meissner currents. However, the experiment of Ref. 6 on SFS junctions shows that the observed Fraunhofcr 
pattern corresponds to < M >= in the F layer. This behavior according to the authors of Ref. 6 may be attributed 
to a multi-domain structure. 

It would be interesting to carry out experiments on S/F structures with non-collincar magnetization in order to 
observe this new type of superconductivity. As follows from a semiquantitative analysis, the best conditions to observe 
the Josephson critical current caused by the TC are high interface transparency (small 7b) and low temperatures. 
These conditions are a bit beyond our quantitative study. Nevertheless, all qualitative features predicted here (angle 
dependence, etc) should remain in a general case when one has to deal with the non-linear Usadel equation. 

Another type of experiments that may detect the triplet condensate is measuring the local density of states. As we 
have shown in the second section, the long-range TC may be detected by measuring the local DoS of the F layers. 

We would like to thank SFB 491 and the German-Israeli Foundation (GIF) for a financial support. 

Note added: After the submission of this manuscript a paper 39 appeared in which a detailed study of the critical 
temperature in FSF structure with noncollincar magnetizations in the F layers has been presented. 
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